A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11(02. 18. 8;02.19.2;03.13.4;11.01.2;11. 19.1) 



ASTRONOMY 
AND 
ASTROPHYSICS 
5.2.2008 



ON 



Anisotropic illumination of AGN's accretion disk by a non 
thermal source : I General theory and application to the 
Newtonian geometry 

G. Henri and P.O. Petrucci 

Laboratoire d'Astropliysique, Observatoire de Grenoble, B.P 53X, F38041 Grenoble Gedex, France 
Received ??; accepted ?? 



(N 



> 
m 
cn 
(N 

o 

ON 

Oh; 

6 



Abstract. We present a model of accretion disk where 
the disk luminosity is entirely due to the reprocessing of 
hard radiation impinging on the disk. The hard radiation 
itself is emitted by a hot point source above the disk, that 
could be physically realized by a strong shock terminating 
an aborted jet. This hot source contains ultrarelativistic 
leptons scattering the disk soft photons by Inverse Comp- 
ton (IC) process. Using a simple formula to describe the 
IC process in an anisotropic photon field, we derive a self- 
consistent angular distribution of soft and hard radiation 
in the Newtonian geometry. The radial profile of the disk 
effective temperature is also univoquely determined. The 
high energy spectrum can be calculated for a given lep- 
ton distribution. This offers an alternative picture to the 
standard accretion disk emission law. We discuss the ap- 
plication of this model to Active Galactic Nuclei, either 
for reproducing individual spectra, or for predicting new 
scaling laws that fit better the observed statistical prop- 
erties. 

Key words: Galaxies: active - Galaxies: Seyfert - Ac- 
cretion, accretion disks - Ultraviolet: galaxies - X-rays: 
galaxies - Radiation mechanisms: non-thermal ~ Scatter- 
ing 



1. Introduction 

With the development of high energy telescopes, it has 
been recognized that Active Galactic Nuclei are the most 
powerful emitters of high energy radiation in the Universe. 
However, the detailed production mechanism is still a mat- 
ter of debate. For radio loud AGN, the detection of very 
high energy radiation, in the GeV (von Montigny et al. 



et al. 1984 ; Dermer & Schlickeiser 1992 ). No such con- 
clusion can be drawn up to now for radio-quiet objects 
such as Seyfert galaxies, since their high energy spectrum 
is apparently cut-off above a few hundred keV (Jourdain 
et al. |1992| ; Maisack et al. |1993| ; Dermer & Gehrels [19951 ). 
Although this radiation could be produced by direct syn- 
chrotron mechanism, it is more often assumed that it 
comes from the Comptonization of soft photons by high 
energy electrons or pairs. Two classes of models have been 
proposed so far: Comptonization by a thermal, mildly rel- 
ativistic, plasma, resulting in a lot of scattering events as- 
sociated with small energy changes, or Inverse Compton 
(IC) process by one or few scattering events from a highly 
relativistic, non thermal particles distribution, which can 
result from a pair cascade. 

Detailed observations in the X-ray range by the Ginga 
satellite have shown that a simple power law is unable 
to fit the X-ray spectrum of Seyfert galaxies. Rather, the 
spectra are better reproduced by a complex superposition 
of a primary power law, with an index a ~ 0.9 — 1.0, a 
reflected component from a cold thick g fluorescent 
Fe Ka line and an a bsorption edge by a warm absorber 
(Pounds et al. |l990| ; Nandra & Pounds |l994[ ). The sec- 
ond and third components could be produced by the re- 
flection of primary hard radiation on an accretion disk 
surrounding the putative massive black hole powering the 



AGN (Lightman & White |l988t George & Fabian |1991 
Matt, Perola & Piro 1991). This has led to consider vari- 
ous geometries where the hot source is located above the 
disk and reilluminates it, producing the observed reflec- 
tion features. The hot source can be a non-thermal plasma 
(Zdziarski et al. 1990), or a thermal hot corona covering 



the disk (Haardt & Maraschi 1991, 1993; Field & Rogers 



1993) 



19951 ) and even TeV (Punch et al. |l992| Quinn et al. |1996D 
ranges, proves the existence of ultrarelativistic particles, 
probably associated with a relativistic jet (e.g. Begelman 



In another context, some observational facts have mo- 
tivated the development of so-called reillumination mod- 
els, where high energy radiation reflected on a cold sur- 
face (presumably again the surface of an accretion disk), 
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produces a fair fraction of thermal UV-optical radiation. 
Firstly, long term observations have shown that for some 
Seyfert galaxies, such as NGC 4151 (Perola et al. 1986) 



and NGC 5548 (Clavel et al. 1992), UV and optical lumi 



nosities were varying simultaneously, and correlated with 
X-ray variability on time scales of months, whereas the 
rapid, short-scale X-ray variability was not seen in optical- 
UV range. This is in contradiction with the predictions of 
a standard, Shakura -Suny aev (SS) accretion disk model 
(Shakura & Sunayev 1973 ), where any perturbation caus- 
ing optical variability should cross the disk at most at 
the sound velocity, producing a much larger lag between 
optical and UV than what is actually observed. Rather, 
these observations support the idea that optical- UV radia- 
tion is largely produced by reprocessing of X-rays emitted 
by a small hot source, the UV and optical radiation be- 
ing emitted at larger distances. The main problem is that 
the apparent X-ray luminosity is usually much lower than 
the optical-UV continuum contained in the Blue Bump, 
whereas one would expect about the same intensity in 
both components if half of the primary hard radiation is 
emitted directly towards the observer and the other half 
is reprocessed by the disk. 

In many cases also, the equivalent width of the Fe Ka 
line requires more impinging radiation than what is actu- 
ally observed if explained by the reflection model (Weaver 
et al. 1995, Nand ra et al. 1997). As an explanation, Ghis- 



ellini et al. (1991), hereafter G91, have proposed that the 
anisotropy of soft radiation could lead to an anisotropic IC 
emission, with much more radiation being scattered back- 
ward than forward. Due to the complexity of their calcu- 
lation, they have restricted themselves to the emission by 
a hemispheric bowl (equivalent to an infinite plane), that 
could model a flared accretion disk with a constant tem- 
perature. The thermal disk-corona model faces the same 
kind of difficulties, for it predicts nearly the same lumi- 
nosity in X-ray and UV ranges. A possible solution could 
imply a patchy corona, a fair part of the UV luminosity 
being emitted by internal dissipation in the disk (Haardt 
& al. 1995). In the cases where X and UV luminosities are 
comparable however, it is difficult to explain very rapid 
X-ray variability as the corona must cover a large part of 
the disk. 

Although the spectral break observed by OSSE around 
100 keV seems to favor thermal models and disprove the 
simplest pair cascade models, such a break could also be 
obtained by a relativistic particles distribution with an 
appropriate upper energy cut-off, such can be provided 
for example by pair reacceleration to avoid pair run-away 
(Done et al. Henri & Pelleticr |l99l|). The aim of 

this paper is to reconsider the reillumination by a non 
thermal, optically thin IC source, taking properly into ac- 
count the disk geometry and the anisotropic distribution 
of photons. We first establish a simple expression to evalu- 
ate the power emitted by a single particle scattering pho- 
tons by Compton mechanism in the Thomson regime in 



an arbitrary soft photon field. The formulae require only 
the computation of the components of the relativistic ra- 
diation tensor, or equivalently the Eddington parameters 
for an axisymmetric field. We then develop a self consis- 
tent model where the emission of the disk is entirely due 
to the reprocessing of hard radiation, produced itself by 
IC process in a hot point source located above the disk. In 
this case a unique angular distribution of hard radiation 
and a unique (properly scaled) disk temperature radial 
profile are predicted. We discuss then the possible phys- 
ical mechanisms for such a situation and its implication 
for the overall characteristics of AGNs, both for individ- 
ual spectra and for statistical properties. We derive new 
scaling laws for luminosity and central temperature as a 
function of the mass. We show that the predictions of the 
model are sensitively different from the standard ones, and 
that they could better explain the observations. 

2. Anisotropic Inverse Compton process 

2.1. Total power emitted by a single particle 

We first establish useful formulae to compute the Inverse 
Compton (IC) emissivity of a particle in an arbitrary pho- 
ton field, in the Thomson regime. We consider the case 
of a relativistic charged particle with mass m, velocity 
V = /3cfco, and Lorentz factor 7 = (1 — /3^)~^/^, in a soft 
photon field characterized by the specific intensity distri- 
bution /i/(fc) . k and feg are respectively the unit vectors 
along the photon and the particle velocity. We assume that 
the Thomson approximation is valid, that is £7 <gC 1 where 
e — hv/mc^ is the soft photon energy in units mc^. In this 
limit, the rate of energy transferred from the particle to 
the photons is: 



dE 



dt 
where 



P+^cTT J h{k) (1 - (3ko.k)dndv 



and 



P_ = CTTT^ j h{k) (1 - pko.kfdQ.dv 



(1) 



(2) 



(3) 



are respectively the power brought by the incident pho- 
tons and carried out by the scattered ones. Here <tt = 
6.65 lO^'^^cm^ is the usual Thomson cross section. 

To transform these expressions, it is useful to consider 

the decomposition of the intensity field /(fc) = j I,j{k) dv 

over the spherical harmonics basis: 



ra- 



I{k) = I h{k) dv=Y, Cl,nYrik) 

-n=-l 
1^0 



(4) 
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where, due to the orthonormahty condition 



the coefficients c/m are given by: 



J J I^{k) Yi""* [k)dVLdv 



(5) 



(6) 



(here * denotes the complex conjugate and 5w the usual 
Kronecker symbol equal to 1 if / = and or else). Note 
that because Iu{k) is real, one has the conjugation rela- 
tionship = C;*^. 

Now one can write k^.k — cos a, where a is the angle be- 
tween the particle velocity and the incident photon, and 
use the following expansion formulae: 



cosa = ^ ^ rr(feo)>r*(fc) 



m— + 1 



cos^ a 



3 

87r 
15 



(7) 



m — — 1 



5] Yr{ko)Yr{k) + ^Y,'{ko)Yo'*{kX8) 



Inserting Eq. (|), (7) and (8) in Eq. and using the 

relation (|5|), one gets finally: 



3 

m— — 1 



4^aT7'[coo(l + y)>? 

m=+l „2 m=+2 



_ = _47raT7'/3[cooYi^o° 



1 -f '"^1 



2/3 

1— — 1 m— — 2 



(10) 



(11) 



Thus the computation of the power emitted in any di- 
rection requires the computation of the 9 components c;™ 
\l = 0, 1, 2; — Z < m < +/} of the radiation field (related to 
the 9 independent components of the relativistic radiation 
tensor). These formula can further be simplified in the im- 
portant case of an axisymmetric field. There c/m = ci6om, 
and the relevant spherical harmonic functions are given 
by: 



Yo'iko) 
Y^ko) 



1 



An 
3 



An 
5 



(12) 
(13) 
(14) 



where fi = cos 6*0 = kg.z, z being the unit vector of the 
vertical axis. Using the Eddington parameters 



H 



I„{k) dfidv ~ 
Iv{k) fidfidif 



1 



An 
1 



Coo 



= Cl0 



(15) 



^-2 



1 



367r^ 



2 . 



one gets finally: 



P+ = AnariJ - l3Hfi] 
P_ = 2naTl^[2J + (3^{J ^ K) 
-AI3H^i + I3^{m - J)ii^] 
dE ^ , , 



dt 



K) 



-2(1 + f3^)Hfi + /3{3K - J)fi'^] 



(16) 
(17) 

(18) 



These expressions appear like simple polynomials of order 
2 in /i, involving only the calculation of the three Edding- 
ton parameters. In the case of a ultrarelativistic particle 
7 1, they take the form: 



P+ 
dE 
'~dt 



2naTl'^JFr,^(p) 



(19) 
(20) 



(9) where we introduce the following notations: 



H 

^ 7 

K 

F^^{^l) = [(3-x)-47?/i+(3x-l)/i' 



(21) 
(22) 



2.2. Emitted spectrum 



Although the total emitted power can be cast into the 
above relatively simple forms, there is no such simplifi- 
cation for the spectrum of the emitted radiation. This is 
because photons with a given energy can be produced by 
different combinations of initial energy, incident angle, and 
scattering angles and thus the exact spectrum depends on 
the detailed form of the soft photon distribution and not 
only on the Eddington moments. A complete calculation 
requires the integration of the Klein-Nishina cross-section 
over photon energies, particle energies and relative angles. 
However, in the case of IC scattering of a single particle 
on a monoenergetic, isotropic soft photon distribution, a 
convenient approximation is often to take a (5-function 



ns5{e' - (e')) 



(23) 



where (e') — —7^65 is the mean energy of Comptonizcd 
3 

photons. 
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In the case of a relativistic distribution, this approxima- 
tion is reasonable if the width of soft photon energy spec- 
trum is much less than the width of the particle energy 
distribution. As we shall see, the soft photon spectrum 
predicted by the present model is close to a blackbody, 
and we will keep this kind of approximation. One can eas- 
ily generalize expression ( p3|) to an arbitrary soft photon 
field by taking the appropriate expression for the Comp- 
tonized photons mean energy. It is obtained by dividing 
the emitted power (Eq. (10)) by the rate of photon scat- 
tering. The latter is given by 



o-T J Iu{k) {hv)-^{l- p.k)dndu 



(24) 



A calculation quite similar to that of the previous para- 
graph gives: 



m=+l 



(25) 



7n— — l 



where dim = 



J J I,,{k) {hv)-^Yi^*{k)dVtdv is calculated 

with the photon number flux instead of the energy flux. 
In the case of an axisymmetric photon fleld again, one can 
simplify this expression using the photon number Edding- 
ton parameters: 



(26) 



J - \jh{k) {h:^)-^dfidi^ ^ iA7:)-^/^dao 
H I^{k) {hv)-^ndiMdv = (127r)-i/2dio (27) 

(28) 



H 



J 



to get: 

n, = AnariJ - f3Hfi). (29) 
The mean photon energy of the emitted radiation is thus: 

Us 

^ [2J + P^{J-K) -4I3HJ1 + P\3K - J)^i^] 
^ 2{J-l3Hn) 

For ultrarclativistic particles, these expressions become 



Ug = 47r(TT'/(l — '7m) 



(31) 
(32) 



J 



where (e^) = = is the mean energy of incident soft photons 
and 



2(1 -W) 



(33) 



is an angle-dependent numerical factor. For an isotropic 
photon distribution, rj =7] = Q and % = 1/3 and one gets 



the familiar result A{fi) = 4/3. Just as in the isotropic 
case, one can approximate the spectrum by a Dirac dis- 
tribution of Eq. (p3|), if most of the emitted energy comes 
from a restricted range of soft photons energy and direc- 
tion. One can expect this to be a good approximation if 
the particle energy distribution is broad enough, so that 
the intrinsic broadening due to photon energy distribution 
is negligible, except near the spectrum energy cut-offs. 

2.3. Emission by a relativistic particles distribution 

The previous formulae can be applied to the case of a 
relativistic particles distribution. For sake of simplicity, 
we will restrict ourselves to the case of an axisymmetric 
distribution /(7,^), which represents the particle number 
(integrated over the volume) per energy and angle cosine 
interval. Axisymmetry is automatically insured at first ap- 
proximation by the cyclotron precession around a small 
magnetic field aligned with the symmetry axis of the radi- 
ation field. For an isotropic distribution, /(7, /i) = n{j)/2, 
where n{"f) is the particle energy distribution. The plasma 
is assumed to be optically thin, such that every particle 
experiences the same radiation field. This point will be 
further discussed in Section (4.4). We assume further that 
the low energy cut-off is high enough to make Eq. (19) - 
(20) valid, whereas the high-energy cut-off is still in the 
Thomson regime. 



2.3.1. The integrated power 

The integrated plasma emissivity can be written 



(34) 



Inserting Eq. (20) yields 

'^^^tJ j F^x(m)/(7,m)7^c?7- 



(35) 



Defining the normalized angular distribution function 

1 n 



7V(72) 



/(7,yu)7 (i7 



where N is the total relativistic particle number and 

(7') J^\fiJ^^^''fh,fih'd^ 



(36) 



(37) 



is the mean quadratic Lorentz factor, one can rewrite this 
expression under the form: 



dP 

dn 



- aTN{j^)Jg{fi)E^^{fi) 



(38) 



The anisotropy of emitted radiation appears thus sim- 
ply as the product of an anisotropy factor of the parti- 
cles distribution g{fJ.) times the anisotropy factor f^^(/i)of 
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the radiation field. For an isotropic particles distribution, 
= 1/2. An interesting case is that of a plasma moving 
relativistically with a bulk velocity /S^and a corresponding 
Lorentz factor 7^, such can exist in super luminal radio- 



sources (Marcowith et al. 1993). The particle Lorentz fac- 



tors in the observer frame 7 and in the plasma rest frame 
7' are linked by the relation: 

7 = r7' (39) 

where T = [7f,(l — PbfJ-)]^^ is the usual Doppler factor. 
Using the Lorentz invariance of /(7,^)/7^, one gets : 

p5 



(40) 



Alternatively, one can express the results in terms of 
mean quadratic Lorentz factor in the plasma frame 
(7'^) = (7')76'(l+/56')toget: 



rIP 

— = aTiV(7 ')r^F,^(M)- 



(41) 



It should be stressed that the above formula holds for the 
emissivity at a given point at rest with respect to the 
observer, such as a small volume of a continuous jet. If 
one follows the plasma in its motion, as in the case of a 
traveling blob, the reception time interval dtn is related to 
the rest time dt by dt = Tdtfi and an extra Doppler factor 
r must be added. In any direction, the emitted power will 
vary as F^. 

2.3.2. Emitted spectrum 

One can also derive the emitted spectrum in a given di- 
rection by integrating the single particle emissivity over 
the particle distribution /(7,/i), that is 



dP 



1 /■^— , dn 



dD.de' 27r 



dtde 



if ill lj)dl 



(42) 



Using the 5 approximation of Eq. (|2^) to integrate over 7, 
and the expressions (31) - (32) , one gets: 



dP 



dQ.de 



- = aT{J-H^Ji){{e,)A{^Ji)r"\'^'^ 



/ 7 = 



(es)A(A*) 



1/2N 



2.4^. Application to a conical photon field 



In this section, we illustrate the simplification brought by 
the present formalism, by redressing the problem of IC 
process in a semi- isotropic photon field. This question 
has already been treated in G91, but the use of non ana- 
lytical integrals requires numerical computation. Here we 
show that their results can be very easily recovered in 



the present formalism, leading to simple analytical expres- 
sions. 

We consider in fact the more general case of a radiation 
emitted in a cone with opening angle 9, with constant spe- 
cific intensity inside the cone and a null intensity outside. 
This would be relevant for example for an isothermal disk 
with a finite radius, or at the vicinity of a star without 
limb darkening. The case = 7r/2 corresponds of course 
to the case studied in G91. In this case, elementary inte- 
gration gives the following Eddington parameters 



J 



cUi: 



-(1-C0S6I) 



J, 



H = -(1 + cose*) 
K = ^(1 + cos 61 + cos^ 6*) 



(44) 



where Uiso denotes the energy density emitted by the full 
sphere as in G91. Equations (19), (20) read now: 



P+ = archil -cos9)[l-^f,] 



(45) 



P_ = aTC^^j'^{l-cos0)[2 + ^{2 + cose + cos^ 



4 ' ""^'^ ' 3 
-2/3^(H-cos6i) +/3V^cos6i(l -I-C0S6')] (46) 

^ ^ _^yc^72/3[^(8-cos0-cos2(?) 
dt 4 3 

-^(1 -|-cos6')(l +/32) +/3/i2cos6i(l -I-C0S6')]. (47) 

The integrals in Eq. (5) of G91, corresponding to cos 9^0, 
appear thus like simple cosine functions. The ratio be- 
tween the (minimal) power emitted in the forward direc- 
tion and the (maximal) power emitted in the backward 
direction is: 

P-{ii = -l) _ 3 + 3/3(l + cose')+/32(l + cos6'-hcos2 



P-{fi = +l) 3-3/3(1-1-^ 



-1-/32(1 + cos6' -f cos^ 



-(48) 



In the limit (3 ~ 1,9 — 7r/2, one gets a factor 7 as found nu- 
merically by G91. One can also easily evaluate the ratio R 
between the total power emitted in the lower hemisphere 
and the upper one by an isotropic monoenergetic particles 
distribution: 



(43) R = 



_ J°^P-dfi _ (2/32 + 3/3(1 -hcos0)+ 6) 
P-d^i ~ (2/32 - 3/3(1 -I- cos9) + 6) ' 



(49) 



with R 



11 + 3cos( 



when /3 — > 1. Once again, one can 



5 — 3 cos 9 

recover G91's result by setting cos 9 — 0, giving R — 2.2. 
The maximally anisotropic case corresponds to cos 6* — 1 
(point source) and gives i? = 7. It is thus the absolute 
maximal ratio between the power radiated in two hemi- 
spheres by an isotropic particles distribution in any pho- 
ton field. 
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I relativistic leptons source 




Fig. 1. The general pictures of the model. We have drawn straight trajectory of a beam of photons, emitted by the hot source 
in a solid angle dQ., and absorbed by a surface ring of the disk dS. We use the notation = cos 6*0 



3. The self-consistent reilluminated disk:the New- 
tonian case 

3.1. Assumptions oj the model 

We consider now a self-consistent model where Inverse 
Compton process takes place on soft photons from the 
accretion disk, which are themselves emitted as thermal 
radiation due to the heating of the disk by hard radiation. 
The disk is modelized by an infinite slab radiating isotrop- 
ically like a black-body at the same equilibrium tempera- 
ture (Fig. |l|). The high energy source is assumed to be an 
optically thin plasma of highly relativistic leptons, at rest 
at a given distance Zq above the disk axis. Its size is small 
enough to be considered as a point source. We consider a 
Euclidean geometry so that there is no curvature of pho- 
ton geodesies. The general relativistic case will be studied 
in an accompanying paper (Petrucci & Henri 1997). The 
particle distribution is assumed to be isotropic. As long as 
the disk emission is concerned, there is no need to specify 
the energy distribution giving rise to IC process. However, 
the spectrum of high energy radiation will be determined 
by this distribution. This will be developed in section ^ 

3.2. The self consistent solution 

We shall see now that the above problem admits in fact a 
unique self-consistent solution in conveniently scaled vari- 
ables. For a given disk emissivity, the power emitted per 
unit solid angle by the high energy source is given by Eq. 
(p8[), where g{^) — 1/2 and the Eddington parameters are 
to be calculated with the disk emissivity. Integrating over 
solid angle, one gets the total high energy luminosity: 



Thus, the power per unit angle can be expressed as: 



dP 

dn 



(51) 



The Eddington coefhcients can be at turn calculated if 
one knows the disk emissivity. Under the hypothesis that 
the disk reprocesses the whole radiation impinging on it, it 



is determined by equating the power absorbed and emitted 
by a surface element dS of the disk at a radius r (cf. Fig. 
{§)■■ 



dP 

F(r)dS = —dn 
dil 



^dPdS 

dnz^ 



(52) 



where ii' = —fi = Zo{r^ + Zq)^^/^ is the cosine of the 
impinging angle of radiation (/i' varies from 1 to 0). 
From equation (pTl) and (&2), one gets: 



Fir) = 



^r,x(-Ai')M 



/3 



(53) 



Under the assumption of isotropic emissivity of the disk, 
one gets the specific intensity due to the reemission by the 
disk at the radius r toward the source: 



F{r) 



where we define the dimensionless parameter 



327r2Z2j' 



(54) 



(55) 



Inserting this expression in the definition of Eddington 
parameters, one gets the following linear system 



1 = C(Go(l -x)+ 4Gi7? + G2(3x - 1)) 
ri = e(Gi(l-x)+4G27? + G3(3x-l)) 
X = e(G2(l - x) + 4G3r/ + G4(3x - 1)) 



(56) 



where 



(50) G„ = 



2(?i + 4) 



n+4 

'max 



m::+„')- (57) 



This homogeneous system admits a non trivial solution 
only if the determinant is set to zero, which gives a cubic 
equation in ^. One can then determine the angular param- 
eters rj and x, and the emitted flux with Eq. (p3|). In the 

case of an infinite slab, one gets simply G„ = — -. 

The numerical solutions of the cubic equation are then 
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<^ 1.449, 48.136 and 7786.45. Only the first one is com- 
patible with the physical constraints ?7 < 1, X < 1. The 
solutions of the system are approximately rj = 0.82288 
and X = 0.69957. 

With these values, the high energy emissivity has the fol- 
lowing universal angle dependence: 



dP 

— in) = it(0.0686 - 0.0982/^ + 0.0328/i2), 
ai2 



(58) 



and thus the illuminated disk has the corresponding emis- 
sivity law: 

Fir) = aT^jfir) 
LtZa 



(Z2 +r2)3/2 
( 

0.0686 + — 



0.0982 



1/2 



0.0328 



(59) 



The total disk luminosity is 

/•OO 

Ldiak = / 2TrrF{r)dr = Lt 

JQ 

~ 0.8086L,. 



377 



(60) 



one can write the emitted spectrum in the universal form: 

T{^') = a^''V4f,^(-^')'/' (65) 
4577^ 



167r4 



r(M') 



1 



dn'.{66) 



It is interesting to compare this spectrum with that of 
a standard SS accretion disk. The flux emitted by a SS 
disk with a null torque at the inner radius rt is given by 



Fir) ^aT^ff.sir) 



2'Kr^ 



r- \ 1/2 
r 



(67) 



For large z/, emitted at large distances, the spectra have 
the same slope, with the dependence (x v^^^. This is 
because the release of gravitational energy and the illumi- 
nation by a central source give rise to the same dependence 
of the energy flux at large distance Fir) oc r~^. 



Inspection of formulae (|5^), (61) and (|67|) show that the 
low frequency luminosity will be the same for both models 
if 



Ldisk _ 2(4 + 37?) f u 

,^0 



-{SS) 



''disk 



(3-x) 



(68) 



For a Schwarzschild black hole with a mass M, ri 
QGM/c^. 



It represents thus the main part of the total bolometric 
luminosity Lt- 

3.3. Disk emission spectrum 

Under the assumption of isotropic blackbody emission, 
without limb darkening, the emitted spectrum, integrated 
over angles, can be written as 



nB^iTef f)2TTrdr 



(61) 



where the effective temperature Te// is defined in equa- 



tion (59) and B^iT) 



2/11/3 



exp 



hi 



ksTeffir) 



- 1 



is the usual Planck source function. Introducing the char- 
acteristic variables: 



/ 3it 



1/4 



\32TTaZ^ 
Vc = ksTc/h 

and defining the following reduced variables: 



(62) 
(63) 



T 



T/Tc 



(64) 



Lrr — — -Li/, 



0.1 



0.01 




1 10 
Reduced frequency v/v^ 



100 



Fig. 2. Spectra of the present illumination model for an infinite 
slab (plain line) and a standard SS accretion disk in reduced 
units (see text for definitions). The 2 spectra are normalized 
in order to have the same low frequency flux 



For comparison, we have chosen Zq = lOGM/c^. Figure 
(H) represents the exact spectrum emitted by the disk in 
our illumination model, with = and fi'max = 1j 

and in the standard SS model. Figure (|^) represents the 
corresponding radial temperature profile for both models. 
As one can see, in the illumination model, the tempera- 
ture stops rising at a distance r ~ Zq whereas it keeps 
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100 F 



10 ^ 




10 100 

radius R/M 




0.1 r 



0.01 



Fig. 3. Radial temperature profile of the present illumination 
model for an infinite slab (plain line) and a standard SS ac- 
cretion disk (dashed line). The 2 curves are also normalized in 
order to have the same asymptotic behaviors for r oo 



Fig. 4. Ratio of high energy to soft energy luminosity as a 
function of inclination angle, in logarithmic scale. The disk is 
modeled by an infinite slab 



growing in the standard accretion disk model, where most 
of accretion energy is released at the smallest radii. As a 
consequence, for the same low frequency flux, the illumi- 
nation model predicts a lower bolometric luminosity, and 
peaks at a lower frequency than the standard SS model. 

3.4- Ratio of high energy to disk emission 

The model predicts also a definite ratio between the high 
energy IC luminosity and the disk thermal emission; the 
total (non intercepted) IC luminosity is emitted in the 
upper hemisphere (0 < /x < 1), whereas the disk luminos- 
ity is equal to the IC luminosity in the lower hemisphere 
(— 1 < /i < 0). From equation (^l|), one gets 



R 



^disk 

1-377/4 

r 



(69) 



One has R ^ 0.237 for an infinite plane. One can also 
evaluate the ratio between the apparent disk and high 
energy luminosities in a given direction of observation /xq- 
It is given by: 



Ri^io) 



n{dP/dn)ic 



disk 



3-Frjx(A^o) 

4//o(4 + 3??)' 



(70) 



This function is displayed in Fig. It presents a 
minimum on the disk axis (/io = 1), both because the disk 
emission is maximal at small inclination angle, and IC is 
minimal because of the soft photon anisotropy. Thus, the 



model predicts that the observed X/UV ratio is smaller 
than 1 for /iQ > 0.2, and can be as low as 0.012. 



4. Application to AGN 

4-.1. The high energy spectrum 

As reminded in the introduction, the high energy emis- 
sion of Seyfert galaxies can be well reproduced by a pri- 
mary power-law spectrum with a spectral index a ~ 0.9, 
superimposed on more complex structures, that can be 
produced by the reflection of the hard X-rays on a cold 
surface. Although the precise modeling of such a reflec- 
tion component is beyond the scope of this paper, it is 
clearly compatible with the above picture. The primary 
power-law emission should be associated with the emis- 
sion from the hot source, and the UV-optical component 
(Blue Bump) associated with the reprocessed radiation, 
together with a Compton backscattered component in X- 
rays (not taken into account in the present model). No- 
ticeably, this power law is exponentially cut-off above a 
characteristic energy of about 100 keV, with some un- 
certainty its precise value. Contrarily to thermal models, 
where this cut-off is related to the temperature of the hot 
comptonizing plasma, we propose to interpret it as a high 
energy cut-off of the relativistic energy distribution. If the 
UV bump is located around 10 eV, this gives an cut-off 
Lorentz factor 70 ~ 10^. Although a detailed model of the 
high energy source is again out of the scope of this work, 
one can note that a model associating pair production and 
pair reacceleration can provide such upper cut-off, to avoid 
catastrophic run-away pair production (Done ct al 
; Henri & Pelletier |^ ) . 



1990 



To account for the high energy cut-off, we will assume that 
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the particles (electrons or positrons) energy distribution 
function (integrated over volume) has the form: 



n(7) = 7Vo7" 



exp 

70 



7r) 



< 7 < 7n 



(71) 



-2 



Inserting this function in Eq. (43), noting that /(7,/i) = 
n(7)/2, one gets the expression of high energy specific 
power: 



dP 
dQ,dv 




(s-3)/2 



-is-l)/2 
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(72) 



2 4 

Reduced frequency v/v^ 



where vq = jQ{es)/h is a high energy frequency cut-off. 

Integrating Eq. (|8|) over all angles, one can also write the Fig. 5. Broad band spectra for various inclination angles, in re- 
total high energy luminosity as: duced units. The values of parameters are the following: s = 3, 

Vo = 10'' 



where Co is expressed as an incomplete gamma function 



Co 



X ^ex.p{-x)dx = r(3-s;7™„/7o;7r, 



One can also evaluate the mean photon energy 



= = Cl hl^r 

J 



where 



8C(4) 



mi)I!::::n^^rd^, 



(75) 



(76) 



4-2. Validity of the point source approximation 

It is obvious that a point source is a convenient, but unre- 

/ alisticx apnroximation of the real geometry of the source, 
/ garnrnao)[T4) . , . „ . 

since it has a zero cross section and a infinite Thomson 

opacity. Rather, for a given number N of scattering par- 
ticles, one gets a minimal size Rmin for the source being 
optically thin; for a homogeneous sphere, it requires 



(79) 



where C('^) is the Riemann Zeta fimction. One can finally which yields: 
compute the photon number anisotropy parameter 77, ap- 
pearing in the computation of A{ii) (cf. Eq. (|33|)): 



This hot source will sustain a solid angle ~ 
T^{Rmin/ . The total number of particles in the opti- 
cally thin regime can be calculated by Eq. (|55|) and (^, 



N 



(80) 



__ H _ 



(77) 



(IT (7^) 

Using the distribution function given by Eq. ([7l|), one gets 



remembering that the specific intensity emitted by the 
disk is supposed to be the Planck function. For /imin = 
and fjimax = 1, one gets Ci ~ 3.41 and rj ~ 0.788 Using the 
reduced variables defined by Eq. (64), one gets finally the 
following expression for the high energy specific power: 



(7^ 



7o 



2 r(3 - 5;7mm/7o;7max/7o) 
r(l - S;7mi„/7o;7ma£i;/7o) 



(81) 



Combining Eq. (|79|), (gO|) and (|8l|), one obtains a esti- 
mate of the minimal radius of the source: 



dP 
dVLdtP 



37o 



s-3 



647rCo 
F^x(A*)exp 



(CiA(Ai))(^-')/'l7(i-^)/2x 



(78) 



_ J_ / 3gr(l - S; -1m.inho; Imax/lo) 

Zq 7o V 2r(3 - S; -yminho] Imax/lo) 



1/2 



(82) 



1/2 



with vq — vq/vc- Noticeably, the high energy cut-off de- 



Figure (|^) represents the contour plot of Rmin/Zo as a 
function of 70 and s, for 7mm = 1 and jmax = 00. As 
we have mentioned, high energy observations of Seyferts 
seem to favour values of 70 — 10^, and s ~ 3. This corre- 



pends on the inclination angle, the smallest angles giving sponds to R/Zq — 0.5, which means that the hot source 



the lowest cut-off. 



is not really point-like, but moderately extended: it could 
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400 500 



yO 



Fig. 6. Contour plot of the minimal angular size of the source 
Rmin/Zo, as a function of the spectral index s and the cut-off 
Lorentz factor 70. The limits of the particle energy distribution 
are 'jmin = 1 and jmax = 00 



15rg, located 
SOvg. These values seem reasonable for a 



be realized by a region with a radius R 
around Zq 

shock in the inner region of a jet emitted by an accre- 
tion disk. Of course, a correct treatment should take into 
account the finite size of the source, but the calculations 
are much more involved if the particles are off-axis: the 
present theory must be considered as a preliminary one, 
and the application to the extended case will be treated 
in a future work. 

4-3. Application to broad-band spectra. 

Figure (||) represents the overall spectra predicted by 
the model for different inclination angles, ranging from 
^0 = 0° to 60°. For all these angles, the X-ray luminosity is 
apparently smaller than the UV bump. Higher inclination 
angles will probably lead to strong absorption through 
the external parts of the disk, presumably a molecular 
torus: in the unification scheme, they would correspond 
to Seyfert 2 galaxies (but see below). The UV/X ratio de- 
pends markedly on the inclination angle: noticeably, the 
smaller ratio correspond to almost face-on objects. The 
upper energy cut-off is also lightly angle dependent, due 
to the A{^) factor in Eq. (79). Face-on objects have the 
lowest cut-off. 

Direct comparison with observations is difficult at this 
stage, because the model does not include other impor- 
tant components, like the absorption edge, the Compton 
reflection feature and the fluorescent Fe Ka line that are 



observed in many Seyferts. Without detailed calculations, 
it can be expected however that these components are 
more pronounced that those predicted by an isotropic il- 
lumination model, the mean enhancement factor being of 
the order of ~ 4 (Eq. (69)). Ginga and more recent 



ASCA observations (Nandra et al. 1997) have found large 



equivalent widths for the Fe Ka line, with a mean value 
around 230 ± 60 eV, but up to SSO^^iao eV for NGC 4151, 



whereas detailed calculations of George and Fabian ( 1991 ) 
and Matt et al. (1991) predicted a value around 140 eV. 
These results depend however on the ionization state of 
the matter (Matt et al. 1993 ) and the assumed iron abon- 
dance (Reynolds et al. 1995 ). Clearly more work must be 
done to clarify this issue. 

The resolved line profiles seen by ASCA can be well fit- 
ted by an illuminated accretion disk in rotation around a 
black hole, with ~ 50 per cent of the line emission orig- 
inating within ~ 20 r^, and ~ 80 per cent within ~ 100 



(Nandra et al. 1997): this is thus compatible with a hot 



source located at a few ten rg above the disk. 
Although a precise prediction of UV/X fiux is difficult to 
assess, we have used the work of Walter & Fink (1993), 
who compared the X and UV flux in a large sample of 
Seyfert 1 galaxies, to make a rough comparison with our 
model. We assumed that the 1375 A flux is close to the 
maximum of the UV bump, whereas the 2 keV flux is not 
strongly contaminated by the soft excess or the reflection 
component, and is representative of the low energy part 
of the hard X-ray spectrum given by Eq. (72). We thus 
consider the apparent UV/X ratio Ra of these two values. 
In our model, the maximum of the UV bump is found 
numerically to be vLq^.max — 0.419 at 17 ~ 4. With the 
geometrical factor due to disk inclination, one gets thus 



_ dPuv 



0.419 



(83) 



The hard X-ray flux at low energy is given by Eq. (72) 
with V Taking s = 3 for simplicity, one gets 



dildU 
, so that: 

Raifi) ~ 9.05 



647rCo 



(84) 



(85) 



Now for a sample of Ng galaxies oriented randomly be- 
tween fj, — and /i = 1, one expects a uniform distribution 



(86) 



One gets the probability distribution of the UV/X ra- 
tio Ra'. 



dP 
dRa 



dfi 



dRa 



9.05Co(3-x-(3x-1)m2) 



(87) 
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Fig. 7. Comparison of the probability distributi on of apparent 
UV/X ratio Ra compiled by Walter and Fink( |l993| ) and the 
predicted values, with 70 = 10^ 

Figure represents the comparison of observed and 
predicted probability distribution of the apparent UV/X 
ratio Ra: for 70 — 10^. To derive the probability density 
dP/dRa, the observed values have been binned into in- 
tervals of width 5, and statistical error bars have been 
added. The agreement is very satisfactory, which can be 
rather fortuitous in view of the approximations used. 
Although the widespread opinion, based on unification 
models, is that Seyfert 1 galaxies are seen nearly pole-on, 
whereas Seyfert 2 are their edge-on counterparts, there is 
some evidence that the reality may be more complex. In 
a statistical study of the 48 Se yfert galaxies from the CfA 
catalog, Edelson et al. ( 1987 ) concluded that all Seyfert 
2, but only one third of Seyfert 1, present an excess at 60 
fim, attributed to thermal emission from dust. Since emis- 
sion from dust is not expected to be highly anisotropic, it 
would imply that the unification model applies only to a 
subclass of objects, which possess an obscuring torus with 
some opening angle. The other subclass could not possess 
this torus, and give only Seyfert 1 galaxies seen under 
any inclination angle. This could explain the presence of 
high inclination angles, i.e. low UV/X ratios galaxies. But 
of course, a large intrinsic or extrinsic absorption could 
change a lot this picture. 

4-4- Scaling laws 

As is well known, a usual assumption for the mass- 
luminosity ratio of AGN is that the bolometric luminosity 
is limited by the radiation pressure to the Eddington lu- 
minosity: 



A-KGMnipC 
ut 



This relationship predicts a linear correlation between 
mass and luminosity for Eddington accreting black-holes. 



L (X AI; the corresponding Eddington temperature is 
given by 



1/4 



cx M' 



^1/4 



(89) 



We show here that the present model predicts different 
scaling laws, if one adds some supplementary assumptions 
on the physics of the hot source. At first sight, all equa- 
tions of the models are linear with respect to the global 
luminosity, so no particular relationship is predicted be- 
tween the luminosity and other parameters like M or Zq- 
However things are different if one considers the micro- 
physics and a realistic geometry of the hot source. First, 
let us consider again the upper energy cut-off of the IC 
spectrum discussed in Paragraph ( ^.4| ) . Observations seem 
to show that all Seyfert galaxies have the cut-off around 
the same value, approximately 100 keV; however, when 
taking into account the reflection component, this esti- 
mate could be somewhat higher, up to 400 keV (Zdziarski 
et al. 1995 ). This is sufficiently close to pair production 



threshold to make plausible the idea that this cut-off is in 
some way fixed by a regulation mechanism to avoid run- 
away pair production. In this case, the cut-off energy is a 
physical quantity, determined by microphysics rather that 
macrophysical quantities. The maximal energy of photons 
produced by IC process is of the order 



^qHi/c ~ constant. 



(90) 



Now it is plausible that the size of the source and its 
distance to the black hole are controlled by the global 
environment responsible for the hot source. Conceivably, 
it could be realized through a strong shock terminating 
an aborted jet. If one makes the (admittedly crude) as- 
sumption that all distances scale like the hole radius r^, 
then one gets R/Zq ^ constant. For a given spectral in- 
dex s, Eq. (|8^ ) predicts 70 ~ constant, and Eq. (|90| ) gives 
i^c ~ constant. Finally one gets the following scaling laws: 



Tc cx i/c oc M 







constant 



(91) 
(92) 



that is , the temperature of the Blue Bump does not de- 
pend on the mass and the luminosity scales like the square 
of the mass. Of course, one could observe substantial vari- 
ations of at least one of these quantities if s or Q varies, 
either by a variation of Zq or a variation of R. Interest- 
ingly, observations seem to corroborate these behaviors: on 
a sample of many quasars and galaxies spanning a large 
range of masses, Walter & Fink (|l99|) found that the 
Blue Bump and soft X-ray excess were approximately in 
the same ratio, although the central masses can differ by a 
factor 10''. This is very hard to explain in the frame of con- 
ventional accretion models. In another study, using a spe- 
cific model to describe the width of Broad Lines from the 
emission by the external part of the disk, CoUin-Souffrin & 
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Joly (1991) have deduced the inclination angles and cen- 
tral masses of a sample of Seyfert 1 galaxies and quasars. 
They found a correlation between mass and luminosity un- 
der the form L cx Af^, with (3 = 1.8 ± 0.6. Ahhough these 
results have been obtained on a limited sample, they are 
compatible with the previous results and clearly different 
from those predicted by the standard accretion disk mod- 
els. A rather intriguing consequence is that there is a max- 
imal mass above which the accretion becomes impossible 
by such a mechanism, where the luminosity predicted by 
Eq. (92) gets higher than the Eddington luminosity. From 
Eq. (62) and (|7^), the total luminosity can be written un- 
der the form: 



2567r6G2 



(93) 



Comparing with Eq. (pq), one finds that Lt ~ L^dd for 



AhC^h'^c nip 2 / y-4 



GAir^GaT 



10 



10 



Zq 
30r, 



V5eV 



M, 







(94) 



(95) 



This value is surprisingly close to that commonly in- 
voked for the most luminous known quasars, which ap- 
parently accrete at a near-Eddington rate. It is however 
unclear how accretion would be stopped for higher central 
masses, since the radiation pressure can be effective only 
if the central engine actually radiates. 

5. Conclusion 

We have shown that a model based on reillumination of 
a disk by an anisotropic IC source could lead to a self- 
consistent picture where the angular distribution of high 
energy radiation and the radial temperature distribution 
of the disk are mutually linked and both determined in 
a single way. The model offers a simple explanation for 
the correlated long term variability of X and UV radia- 
tion, the short term variability of X-rays non correlated 
with UV variations, and the apparent X/UV deficit that 
seems contradictory with simple reillumination models. In 
its simplest form, it predicts a unique shape of disk spec- 
trum and a X/UV ratio depending only on the inclina- 
tion angle. The predicted values are in good agreement 
with observations. In an accompanying paper (Petrucci & 
Henri 1997), we shall study the influence of relativistic 



corrections on the above scheme to account for the gravi- 
tational effect of a Schwarzschild or a Kerr black hole. A 
precise comparison with real spectra should also include 
other components, such as a reflection component and a 
fluorescent Fe Ka line. This is deferred to a future work. 
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